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We show that thermally activated interstitial and vacancy defects can lead to first order melting 
of a vortex lattice. We obtain good agreement with experimentally measured melting curve, latent 
heat, and magnetization jumps for YBa2Cu307_i and Bi2Sr2CaCu2 0g. The shear modulus of the 
vortex liquid is frequency dependent and crosses over from zero at low frequencies to a finite value 
at high frequencies. We also find a small surface tension between the vortex line liquid and the 
vortex lattice. 



I. INTRODUCTION 



It has been experimentally established that below 
a critical value of the magnetic field, vortex lat- 
tices undergo a first order transition in clean high 
temperature superconductors [0J21. This has been 
seen in YBa 2 Cu 3 7 _,5 (YBCO) [§f), Bi 2 Sr 2 CaCu 2 8 
(BSCCO) and (La 1 _ 2 ,Sr 2 ,) 2 Cu0 4 |j§. Evidence 

for first order phase transitions comes from latent heat 
measurements || and peaks in the specific heat 0-|9)as 
well as jumps in the resistivity J3|,[n],[^| and in the magne- 
tization |0Mm,O . This transition is generally accepted 
as a melting transition from a vortex solid to a vortex 
liquid. 

A great deal of theoretical work fil|-|24|] has helped to 
establish that there is a first order melting transition. 
Brezin, Nelson and Thiaville [H showed that including 
fluctuation effects in Abrikosov's mean field theory of the 
flux lattice transition would drive the transition first or- 
der. The advent of the high temperature superconduc- 
tors and the subsequent experimental indications of vor- 
tex lattice melting sparked intense theoretical activity. 
Early analytic efforts used the Lindemann criterion [fl5| , 
though such an approach could not show that the tran- 
sition was first order. Studies of the mapping between 
vortex lines and the world lines of a 2D system of bosons 
have suggested a first order transition from an Abrikosov 
lattice to an entangled vortex liquid P5| , |26U241 . Numeri- 
cal simulations [^6[- p3|j27| have been able to show that the 
melting is first order by calculating quantities such as the 
magnetization jump [O^] and the delta function in the 
specific heat pO| , pT| , p3| T^ However these simulations were 
done in the limit of high magnetic fields <C a Q <C X a b 
where t; a b is the coherence length, a D is the spacing be- 
tween vortices, and A a fc is the penetration depth. Most 
p6| , p7 19 -2^] assumed that the magnetic induction B was 
spatially uniform and thus neglected the wavevector de- 
pendence of the elastic moduli. This has made quanti- 
tative comparison with experimental data difficult, espe- 
cially in the case of BSCCO whose vortex lattice melts 



at low fields. In addition the mechanism for vortex lat- 
tice melting is still not well established. There have been 
suggestions that topological defects |28|] such as vortex 
loops |22|,p3[ , vortex-antivortex pairs |2C|j2l[| , free discli- 
nations |L8| , and dislocations |2Sf| may play a key role in 
triggering melting. 



A. Melting Scenario 

In this paper we show that melting can be induced 
by interstitial and vacancy line defects in the vortex lat- 
tice which soften the shear modulus cqq. This softening 
makes it easier to introduce more defects and increases 
the vibrational free energy. The increased vibrations ul- 
timately lead to melting. There is good agreement with 
the experimental curve of transition temperature versus 
field, latent heat and magnetization jumps for YBCO and 
BSCCO. Using a viscoelastic approach, we show that the 
shear modulus is frequency dependent. At zero frequency 
the vortex liquid cannot sustain a shear while at high fre- 
quency the liquid has a finite shear modulus. Since we 
can calculate the free energy for both the lattice and the 
liquid at melting, we have estimated the surface tension 
between a vortex line liquid and a vortex solid. 

Let us describe our scenario for melting. Our ap- 
proach follows that of Granato |30| as well as previous 
work which showed that defects can lead to a first order 
phase transition J3l]]. We start with a vortex lattice in 
a clean layered superconductor with a magnetic field H 
applied perpendicular to the layers along the c-axis. We 
consider the vortices to be correlated stacks of pancake 
vortices. We will assume that the transition is induced 
by topological defect lines, i.e., vacancies and intersti- 
tials. In a Delaunay triangulation |p2[ a vacancy or an 
interstitial in a triangular lattice is topologically equiv- 
alent to a pair of bound dislocations |lq] as well as to a 
twisted bond defect |35| ]. High temperature decoration 
experiments |33| and Monte Carlo simulations Jl£| have 
found such defects to be thermally excited. The intro- 
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duction of these defects softens the elastic moduli. Since 
the energy to introduce interstitials and vacancies is pro- 
portional to the elastic moduli, softening makes it easier 
to introduce more defects. The softening also increases 
the vibrational entropy of the vortex lattice which leads 
to a melting transition. The transition is driven by the 
increased vibrational entropy of the vortex lines of the 
lattice, and not by the entropy of the wandering of the 
defect lines. In fact Frey, Nelson and Fisher |j4j showed 
that a phase transition driven by the entropy of wan- 
dering flux lines occurs at a much higher magnetic field 
than what is observed experimentally. In the vicinity 
of the experimentally observed first order phase transi- 
tion, wandering in the transverse direction by more than 
a lattice spacing is energetically quite costly and there- 
fore rare. Such flux line bending also makes dislocations 
fS5| , j36f energetically costly at low dislocation densities. 
(The energy scale is set by e Q s 0^]. Here s is the in- 
terplane spacing and e G , the energy per unit length of a 
vortex, is given by e c = (0 o /47rA a fc) 2 where </> Q is the flux 
quantum and X a b is the penetration depth for currents 
in the ab plane. For example, for YBCO e a s ~ 650 K at 
T = 70 K and for BSCCO e Q s - 550 K at T = 60 K @- 
Note that e s > T.) 

The first order transition is nucleated in a small re- 
gion by a local rearrangement of existing line segments. 
Slightly above the melting temperature T m a vortex line 
can distort and make an interstitial and a vacancy line 
segment that locally melt the solid. This is the analog of 
a liquid droplet which nucleates melting of a crystal. The 
role of the surface tension is played by the energy to con- 
nect the interstitial segment to the rest of the vortex line. 
This connection can be a Josephson vortex lying between 
planes or a series of small pancake vortex displacements 
spread over several layers. When the length i of the in- 
terstitial and vacancy segments equals the critical length 
£ c , the energy gained by melting equals the energy cost 
of the connections. When I > £ c , it is energetically fa- 
vorable for the defect segments grow to the length of the 
system. We are ignoring the surface tension associated 
with the surface parallel to the c-axis. We shall show 
later that this is quite small. 



II. FREE ENERGY 

To study melting we assume that we have a vortex 
lattice with interstitial and vacancy lines extending the 
length of the lattice. Our goal is to find the free energy 
density as a function of the concentration n of defect 
lines. The free energy density is 



/ Jo ~t~ fw ~t~ fvib fu 



(1) 



where f a is the free energy density of a perfect lattice, /„, 
is the work needed to introduce a straight interstitial or 



vacancy line into the lattice, f v n, is the vibrational free 
energy density of the system, and f wan is the free energy 
due to the wandering of the defect lines over distances 
large compared to the lattice spacing. We now examine 
these terms in detail. 

/o, the free energy density of a perfect rigid flux lattice, 
is given by the London term p3,B8[: 



fo = 



B 2 
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where B is the spatially averaged magnetic induction, £ a {, 
is the coherence length in the ab plane, and rj is 0.130519 
for a hexagonal lattice and 0.133311 for a square lattice 
|p4j . For B near H C 2, f is given by the Abrikosov free 
energy SI 



(H c2 - Bf 
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(3) 



where the Ginzburg-Landau parameter k = X a b/S, a bi and 
the Abrikosov parameter Pa is 1-16 for a triangular lat- 
tice and 1.18 for a square lattice. 

To calculate f v n,, we follow ref. [^0|. We denote 
the displacement of the z^th vortex pancake in the nth 
plane from its equilibrium position by u(n, r„) where 
u = (u x , u y ) and the pancake position r = (r x , r y ). The 
Fourier transform u(k, q) = ^ n!/ u(n, r„) exp[i(k • r„ + 
qn)]. k = (k x ,k y ) and q is the wavevector along the 
c-axis. f V ib = — (ksT/V) In Z v n> where V is the volume 
and the vibrational partition function Z V ib is given by 



Z v ib — 



-T c ilk B T 



n 

k,q>0,i 



dufl(zk, q)duj(i\i, q) 



(4) 



where we have divided by the area 7r£ 2 b of the normal 
core of a pancake fl0|| . u r and uj are the real and imagi- 
nary parts of u(k, q) and ie{x, y}. The elastic free energy 
functional associated with these distortions is 



1 



(5) 



kg ij 



where i and je{x,y}, the volume per pancake vortex is 
v Q = s<p /B, and s is the interplane spacing. The k 
sum is over a circular Brillouin zone K% = 4wB/<j) . The 
matrix aij is given by 



c B kikj + (c 66 fc 2 + c 44 Q 2 )(5i 



(6) 



where cb, cqq, and C44 are the bulk, shear, and tilt mod- 
uli, respectively, cb — en — cqq for a hexagonal lattice. 
Q 2 = 2(1 — cosgs)/s 2 . Diagonalizing ay leads to 2 eigen- 
values: 



A e (kq) = cnfc 2 
A t (kq) = c 66 k 2 



C44Q 2 
caQ 2 



(7) 
(8) 
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where A is the diagonal matrix, the subscript £ denotes 
longitudinal and t denotes transverse. Using this leads 
to 



(9) 



kg i=t,t 



where ie{£,t}. After integrating over u in the re- 
maining sums over k and q are converted to integrals: 



4tt J_„ /s 2tt \v £ b Ai 



(10) 

where the volume of the sample is set to unity. The 
integrals in (|l0|) are done numerically. At low fields 
(b = B/H C 2 < 0.25), the elastic moduli are given by 
"14T 
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where A c is the penetration depth for currents along the 
c-axis, 7 = Xc/Xab is the anisotropy, and £ = 1. At 
high fields (b > 0.5) c 66 is altered by the fac- 

tor 0.5k" 2 )(1 - 6) 2 (1 - 0.586 + 0.2% 2 ) and 

the penetration depths in en and C44 are replaced by 
A 2 = A 2 /(l — b) where A denotes either A a 6 or A c . 
In addition the last two terms of C44 are replaced by 
B(f> / (16ir 2 A 2 ). These replacements guarantee that the 
elastic moduli vanish at H c2 . For YBCO the temper- 
ature dependence of the penetration depths and coher- 
ence lengths are given by A(T) = A(0)(1 - (T/TJ)^ 1 / 3 
@ and U(T) = U(0)(1 - (T/T,))- 1 / 2 , respectively. 
For BSCCO whose melting field is two orders of mag- 
nitude below H c2 , A 2 (T) = A 2 (0)/(1 - (T/T c ) 4 ) and 
e ab (T) = &(0)/(l - (T/T c ) 4 ) gj. 

The free energy density f w due to the energy cost of 
adding a vacancy or interstitial vortex line is difficult to 
calculate accurately However, we can write down 

a plausible form for f w by noting that a straight line de- 
fect parallel to the c-axis produces both shear and bulk 
(but not tilt) distortions of the vortex lattice. For ex- 
ample, if a defect at the origin produces a displacement 
u that satisfies V • u = v S(r)/s where 6(r) is a two di- 
mensional delta function, then it Q (k) = ik a /k 2 p4j , E5j , 
Inserting this in (|^), we find that f w — (cgg+5e)/2 where 
cb = J2k CB ( ( l = 0, k). Generalizing this to allow for a 



more complicated distortion and for a concentration n of 



line defects, we write 30 



dn(a\CQQ + a 2 c B ) 



(12) 



where ati and a 2 are dimensionless constants. We expect 
the isotropic distortion to be small, i.e., a 2 <C 1, and the 
shear deformation to dominate, i.e., a± ^> a 2 . Integrat- 
ing over n allows the elastic moduli to depend on defect 
concentration. We will assume that cb is independent of 
n since we believe that the bulk modulus of the vortex 
solid is roughly the same as that of the liquid phase. To 
find C66 (ri) |30l , we use its definition 



c 66 = d 2 f/de 2 



(13) 



where e is the shear strain. Assuming that cb has negli- 
gible shear strain dependence, we find 



or 



cm(n) = c 66 (0) + ai (8 c 66 (n)/<9e )dn 
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(14) 



(15) 



If we shear the lattice in the ab plane along rows sep- 
arated by a distance d, the system must be unchanged 
if the displacement is equal to a lattice spacing. This 
is a result of the discrete translational symmetry of the 
lattice. The shear modulus should reflect this discrete 
translational symmetry and therefore must be periodic in 
displacements equal to the lattice constant a = y/ <f> /B. 
We describe this with the simplest even periodic function: 



c&g{u) = cgq(u — 0) cos(2wu/a ) 
= cqq(s = 0) cos(27r cfe/a ) 



(16) 



where the shear strain e = u/d. Notice that this ex- 
pression goes beyond the usual harmonic approximation. 
Then taking the second derivative of eq. (|l6|), we obtain 



d 2 c 66 {n)/de 2 = -(3c 66 (n) 



(17) 



where (3 = An 2 d 2 /a 2 . Combining this with (y_5|), we ob- 
tain 



ceeO*) = c 66 (0) exp(-«i/3n) 



(18) 



where C66(0) is given in eq. (|l l|) . Thus the shear mod- 
ulus softens exponentially with the defect concentration 
n. This softening lowers the energy cost to introduce 
further defects, and increases the vibrational free energy 
f V ib when C66 (n) is used in a^- . Substituting cqq (n) in eq. 
(pTsj) into our expression (|l2|) for f w yields 



fu 



c 66 (n = 0) 



(1 



-cti/3n 



) + a 2 CBn 



(19) 
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The softening of the shear modulus with increasing 
defect concentration is well known in the case of atomic 
lattices }i6|. There it has been shown both experimen- 
tally |47-|49[| and theoretically |5(J that interstitials can 
substantially soften the elastic constants with the largest 
change being in the shear modulus. Linear extrapola- 
tion of the experimentally measured change of the shear 
modulus of copper would imply that the lattice becomes 
unstable for a concentration of about 3% interstitials 
p6[ . An example of how interstitials can soften the 
shear modulus is illustrated in Figure El Here we show 
a triangular lattice where an interstitial forms a dumb- 
ell aligned in the < 010 > direction by sharing a site 
with another atom or flux line. Dumbbell displacements 
along the < 100 > direction introduce a string-like libra- 
tional resonance mode consisting of displacements along 
the < 110 > directions. This mode couples strongly to 
an external shear stress and results in softening of the 
shear modulus |^l| . 

The last term we need to consider is f wan , the free 
energy due to the wandering of the defect lines over dis- 
tances large compared to the lattice spacing. We can 
estimate f wan with the following expression |34) 



k B T 



ln(m^) 



(20) 



where mi = 3 for a triangular lattice (BSCCO) and 
Tji£ = 4 for a square lattice (YBCO). i z can be thought 
of as the distance along the z-axis that it takes [j|2| for 
the defect line to wander a transverse distance of one lat- 
tice spacing a . To go from one vacancy or interstitial 
site to the next, the defect line segment must jump over 
the barrier between the two positions. This gives £ z a 
thermally activated form: £ z ~ l Q exp(— E/ksT), where 
£ Q w a Q (ei/es) 1 ^ 2 and E ^(eies) 1 ^ 2 ■ £i is the line 
tension and is given by ei ~ (eo/7 2 ) hi(a /£ Q ;,). Numeri- 
cal simulations |34|j4^ ] indicate that the barrier height e# 
is small and we use £b = 2.5 • 10 _3 e o . f wan itself is quite 
small compared to the other terms because of the high 
energy cost of vortex displacements. For example, at the 
transition f wan is about two orders of magnitude smaller 
than f w or Thus the transition is not driven by a 

proliferation of wandering defect lines because near the 
transition the high energy cost of vortex displacements is 
not sufficiently offset by the entropy of the meandering 
line [@. 

Before we plot / versus n, we note that the difference 
between B and H is negligible for YBCO but can be a 
significant fraction of the melting field H m for BSCCO. 
To find the value of B to use in the Helmholtz free energy 
density /, we minimize the Gibbs free energy density G, 
i.e., dG/dB = where G = f - B • H/4tt. Since the 
concentration dependence of B is negligible, we find B 
for n — for each value of H and T. Typical plots of 
A/ = f(n) — f(0) — f w + Af v n, versus n are shown in the 



inset of figure 1. The double well structure of A/ is char- 
acteristic of a first order phase transition. The equilib- 
rium transition occurs when both minima have the same 
value of A/. We associate the minimum at n = with 
the vortex solid and the minimum at finite n with the 
vortex liquid. The defect concentration at the transition 
is only a few percent. Eq. Jl8| ) implies that a finite value 
of n yields a finite value for the shear modulus, e.g, for 
n = 5%, c 66 (n) - 0.2c 66 (0) for BSCCO. Previous work 
interpreted this to mean that the lattice did not melt. 
However, they did not appreciate the fact that the shear 
modulus is frequency dependent and the cqq used here 
is the high frequency response. At high frequencies it is 
the elastic response which dominates and this is what en- 
ters into the expression for the free energy. For a liquid 
the low frequency response is dominated by viscosity so 
that the zero frequency shear modulus is zero. We will 
elaborate more on this later. 



III. FITS TO EXPERIMENTAL DATA 

A. Melting Curve 

In Figure ^| we fit the experimental first order transi- 
tion curves in the H — T plane using 2 adustable param- 
eters: ai and ct^. As expected, a\ ^> cti and «2 "C 1 (see 
Figure ||). The geometrical quantity /? can have several 
values for a given lattice structure, depending on which 
planes are sheared. We choose (3 = ir 2 tan 2 (j> where cf> is 
the angle between primitive vectors. Decoration experi- 
ments on BSCCO find a triangular lattice [|3), so we use 
<f> = 60°. For YBCO we choose 4> = 44.1° which is very 
close to a square lattice which has <fi — 45°. Maki |5^] has 
argued that the d-wave symmetry of the order parame- 
ter yields a square vortex lattice tilted by 45° from the 
a— axis. Experiments fSJ [57]] on YBCO find (j> ranging 
from 36° to 45°. 



B. Magnetization and Entropy Jumps 

We can calculate the jump in magnetization AM at 
the transition using AM = —dAG/dH\T=T m where 
AG = G(ni) — G(n — 0). Here rig is the defect concen- 
tration in the liquid at the melting transition. The jump 
in entropy As is given by As = —v dAG/dT\H=H m 
where As is the entropy change per vortex per layer. 
The results are shown in Figure [| We have checked 
that our results satisfy the Clausius-Clapeyron equation 
As = —(v AB/Air)dH m /dT. We obtain good agreement 
with experiment for YBCO and the right order of magni- 
tude for BSCCO. The difference between theory and ex- 
periment in the temperature dependence of the entropy 
and magnetization jumps for BSCCO may be due to the 
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decoupling of the planes ]58|-|62]] . This enhances the ther- 
mal excursions and hence the entropy of the pancake vor- 
tices |i3]-[35| . Decoupling may be brought about by other 
types of defects such as dislocations which we will discuss 
in section 5B. 



C. Lindemann Criterion 

We can compare our results with the Lindemann cri- 
terion by calculating the mean square displacement < 
\u\ 2 > at the transition using eq. (||): < \u\ 2 >= 
~{2k B T/v ) X) Q k q 9 m Z mb /dA(a\iq) where A is given by 
(||) and a labels the 2 eigenvalues. Defining the Lin- 
demann ratio cl by c\ =< \u\ 2 > /a 2 , we find that 
c L 0.25 for YBCO at H m = 5 T and that c L « 0.11 
for BSCCO at H m = 200 G. Here we have used the same 
values of the parameters that were used to fit the phase 
transition curves in Figure ||. These values of cl are con- 
sistent with previous values |l],|2|,[l5| . 

D. Hysteresis 

Experiments have found little, if any, hysteresis 
P, [l0|]l1| . This is consistent with our calculations. We 
can bound the hysteresis by noting the range of temper- 
atures between which the liquid minimum appears and 
the solid minimum disappears. Typical values for the 
width of this temperature range are 300 mK for YBCO 
at H = 5T and 1.3 K for BSCCO at H = 200 G. An- 
other measure of the hysteresis can be found in the plots 
of Af versus n (see inset of Fig. ||). The barrier height 
Vb between the minima is low (Vbv ~ 30 mK) which is 
consistent with minimal hysteresis. 

E. Loss of Superconducting Phase Coherence 



IV. SURFACE TENSION 



The vortex line wandering renormalizes the coupling 
between the planes in the liquid phase, so it is difficult 
to estimate the surface tension parallel to the afe-planes 
which is primarily due to Josephson vortices. However, 
since we have expressions for the free energy in both the 
liquid and solid phases, we can estimate the surface ten- 
sion parallel to the c-axis along the melting curve. We 
imagine a plane interface parallel to the c-axis between 
the vortex liquid and the vortex lattice phases. The sur- 
face tension a is given by 

/oo 
(G(x) - G )dx (21) 
-oo 



where G(x) is the Gibbs free energy as a function of posi- 
tion and the constant G is the Gibbs free energy far away 
from the interface, for example within the solid phase 
where n = 0. (At melting the vortex liquid and solid 
phases coexist because they have the same bulk value for 
the Gibbs free energy.) In the interface region the defect 
concentration n changes from zero in the solid phase to 
a finite value in the liquid phase. Let us assume that in 
this region the concentration gradient dn/dx is a constant 
n /a where n a is the concentration of defects in the bulk 
liquid phase. Here we are assuming that the width of the 
interface is of order a vortex lattice constant a a . Then 

o — —— I (G(n) - G a )dn (22) 
n a Jo 



In going from the normal metallic phase to the vortex 
solid, two symmetries are broken: translational invari- 
ance and gauge symmetry which produces the supercon- 
ducting phase coherence along the magnetic field. In the 
liquid, longitudinal superconductivity is destroyed by the 
wandering and entanglement of the vortex lines. Even 
though line wandering is energetically costly and there- 
fore rare, it does occur. As a result, the correlation length 
along the c-axis will be quite long and of order £ z . This 
is consistent with measurements in YBCO of the c-axis 
resistivity which find that there is loss of vortex veloc- 
ity correlations for samples thicker than 100 /im [|66|-|68|. 
For an infinitely thick sample, the loss of longitudinal 
superconductivity coincides with the melting transition 
p9[ . This agrees with experiments which indicate that 
the loss of superconducting phase coherence along the 
c-axis coincides with the first order transition [B6H6q| . 



This is an integral of the area under the barrier between 
the solid and liquid phases in the plot of the Gibbs free 
energy versus defect concentration (see inset of Fig. ||). 
Using the values for T and B along the melting curve, 
we find the surface tension given in Fig. ^. The de- 
pendence of the surface tension on the melting temper- 
ature T m reflects that of the barrier height Vb on T m . 
The order of magnitude of the surface tension is given 
by a ~ V B /sa a , and as a result, the values are quite 
small. For example, at T=60.24 K and 5=202.28 G, 
cr=0.015 K/sa Q for BSCCO, and at T=80.9184 K and 
ii=6.4807 T, a = 7.26 x lQ- 3 K/sa for YBC O, where 
s is the interplane spacing and a a — yj <j> /B. We be- 
lieve these are correct order of magnitude estimates for 
the surface tension since the small values of the barrier 
height is consistent with the small amount of hysteresis 
found experimentally 0^0,0). 
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V. VORTEX LIQUID 

A. Viscoelastic Behavior 

We now discuss the viscoelastic behavior of the vor- 
tex liquid. As we mentioned earlier, the shear modulus 
is frequency dependent. The low frequency response to 
a shear stress is flow and this is characterized by a vis- 
cosity 77. At high frequencies the response is elastic and 
the crossover between the two occurs over a narrow fre- 
quency range so that the shear modulus as a function of 
frequency is rather like a step function. This behavior 
can be simply modeled using the Maxwell model |7(J for 
viscoelasticity in which a massless spring is damped by a 
viscous force. The rate of shear strain s is given by 



c 66 (lo 



a 
V 



(23) 



where a is the shear stress, a is the time derivative of the 
shear stress, and u> is the frequency. Using the Maxwell 
relation for the relaxation time r = r)/cee(u! = 00) and 
the definition of the frequency dependent shear modulus 
cee(u) = cr(u>)/s(u>), we find that the real part of the 
shear modulus is given by 



c 66 (w) = c 66 (uj = 00) 



1 + w 2 r 2 



(24) 



Notice that cqq(oj = 0) = which confirms that the vor- 
tex liquid cannot sustain a shear stress. At high frequen- 
cies Cee(^) is given by cqq{oj — 00). To estimate the 
crossover frequency we need to estimate the viscosity. 

We are interested in the shear viscosity which arises 
from the interactions between vortices. There are other 
sources of viscosity. For example a single moving vor- 
tex line experiences a viscous drag due to the normal 
electrons in the core which produce resistance when they 
move with the vortex. This is described by the Bardeen- 
Stephen model S|. There is also viscosity which arises 
from pinning; we will ignore this contribution since we 
are considering a clean lattice. We can obtain a simple 
estimate for the viscosity following the approach of Dyre, 
Olsen and Christensen jnj . Shear flow occurs when some 
vortices push past other vortices. The viscosity 77 is given 
by 



r] = t] exp 



AF(T) 
k R T 



(25) 



where the prefactor 770 is the viscosity of a single noninter- 
acting vortex line and is given by the Bardeen-Stephen 
relation [ |38| r\ ~ 4> H C 2/ p n c 2 where p n is the normal 
state resistivity and c is the speed of light. In ( p5| ) AF(T) 
is the activation energy which is identified with the work 
done per vortex pancake in shoving aside the surrounding 
vortices. The elastic energy associated with distorting 



the vortices is given by (||); we identify AF(T) with T e i 
at the maximum distortion u produced by the shoving. 
The actual form of the distortion is difficult to determine 
analytically. We will assume that the dominant contri- 
bution comes from tilt and shear; and that there is no 
change in the density so that the contribution from the 
bulk modulus can be ignored. In BSCCO in the vortex 
liquid, the planes are decoupled and the tilt modulus can 
be ignored. In YBCO in the liquid the correlation along 
the c-axis can be quite long as we discussed earlier. In 
this case the distortion involves various wavevectors; the 
wavevector dependence of the tilt modulus C44 is such 
that at small q, C44 is roughly comparable to the high 
frequency shear modulus cq§(u = 00) and at large q, 
C44 <C cgg(w = 00). (Since we are considering the elastic 
response, it is appropriate to consider the high frequency 
shear modulus.) So as a crude estimate we will assume 
the displacement is pure shear and write |72[| 



AF(T) = c 66 (w = 00, n, T)V C 



(26) 



where V c is the volume change due to shoving and rear- 
ranging vortices. Since AF(T) is the energy per vortex 
pancake, V c is some fraction 8 of the volume v per vor- 
tex pancake, i.e., V c = 8v . c§§(uj = 00, n, T) is given by 
eq. jl8|). In Fig. we show the reduced viscosity 77/77,, 
along the melting line for both YBCO and BSCCO with 
8 = 1. As one can see, interactions enhance the viscos- 
ity 77 over the noninteracting viscosity 77,, by a factor of 
2 or less. This is because AF(T m )/kBT m < 1 along the 
melting line. 

Using our estimate of the viscosity and eq. (pi]), we 
can calculate the frequency dependence of the shear mod- 
ulus cee(uj) as shown in Fig. ^| for a defect concentration 
of 5%. As expected the shear modulus has the shape 
of a step function; it is zero at low frequencies and rises 
quite sharply to its infinite frequency value cee(uj = 00). 
Notice that for BSCCO the crossover frequency is a few 
MHz and for YBCO it is a few GHz. Since 1 K corre- 
sponds to 20 GHz, this means that we made an excellent 
approximation in setting cqq = cq§(uj = 00) in the free 
energy density / in eq. (El). 



B. Dislocations 

Any theory that tries to describe melting has to con- 
tain two main ingredients: a satisfactory description of 
both the solid and the liquid phases and a mechanism by 
which the system goes from one to the other. The diffi- 
culty has always been to describe two phases with wildly 
different properties within the same framework. In the 
present work, we have achieve this by viewing the liquid 
as a solid with a finite concentration of vacancies and in- 
terstitial. Clearly this is an approximation. If one wishes 
to describe a liquid as a solid with defects, other types 
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of defects have to be taken into account as well. This is 
particularly true of dislocations. While we do not believe 
that thermally excited dislocations play the key role in 
mediating the melting transition because of their high en- 
ergy cost, we believe that they will proliferate as soon as 
the vacancies and interstitials appear. There are two rea- 
sons for this. First interstitial vortex lines are usually at- 
tractive p4j and can aggregate to form dislocations that 
extend the entire length of the lattice parallel to the c- 
axis. The same is true for vacancies. In particular, Olive 
and Brandt p4| ] have done numerical simulations on line 
defects which were at least 5 lattice spacings apart. They 
found that both centered and edge interstitials J7^] are 
attractive if \ a b/o, > 1- For X a b/a — 0.25, edge inter- 
stitials were attractive for separations less that 10 lattice 
spacings and repulsive at larger distances while centered 
interstitials were repulsive at distances larger than 5 lat- 
tices spacings. Vacancies were attractive in all cases. 

The second reason is that the substantial softening of 
the shear modulus brought about by the vacancies and 
interstitials reduces the dislocation core energy as well 
as the elastic energy of creating dislocations. For exam- 
ple, the core energy of a z-directed dislocation goes as 
C66& 2 and the core energy of a screw dislocation goes as 
^/cggC44b 2 where b is the magnitude of the Burger's vector 
p6| . In addition the long-range interaction between dis- 
location loops is mediated by the strain field and depends 
on C66 p6| ]- Once the dislocations have proliferated, the 
method and results of Marchetti and Radzihovsky |36| ] 
can be used to provide a more detailed and accurate de- 
scription of the liquid side of the transition. For example 
they show that when dislocations at all length scales are 
present, the shear modulus vanishes in the long wave- 
length limit p5| , p6[ |. Work in this direction is in progress. 

VI. SUMMARY 

To summarize we have presented a model for the melt- 
ing of a vortex lattice into a vortex liquid. The melting 
transition is induced by a few percent of vacancy and in- 
terstitial vortex lines that soften the shear modulus and 
increase the vibrational entropy. The increased vibra- 
tional entropy leads to melting. We obtain good agree- 
ment with the experimentally measured curve of transi- 
tion temperature versus field, latent heat, and jumps in 
magnetization for BSCCO and YBCO. The Lindemann 
ratio c L is - 11% for BSCCO and - 25% for YBCO. 
The hysteresis is small. We find a very small surface ten- 
sion between the vortex solid and the vortex liquid along 
an interface parallel to the c-axis. The shear modulus is 
frequency dependent; it is zero at lu = and plateaus at 
higher frequencies to its infinite frequency value. 

We thank Andy Granato and Lev Bulaevskii for helpful 
discussions. This work was supported in part by ONR 
grant N00014-96-1-0905 and by funds provided by the 
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FIG. 1. Dumbbell interstitialcy configuration in a triangu- 
lar lattice introduces a string-like librational resonance mode 
that couples to external shear stress and softens the shear 
modulus. The atoms or flux lines involved in this mode are 
shown with open circles while the rest of the sites are denoted 
by solid circles. 



8 



Temperature [K] (YBCO) 

80 84 88 



q 500 
O 

00 400 



W 300 
CO 

CD 



CD 



, 200 
100 




It) 

* a o 

+ < -10 



• T<T m 


/t=t; 







0.00 0.02 0.04 0.06 0.0 



> Theory 
- Experiment 



CD 
cn^ 

^< 

CD 

o 
o 



50 70 90 

Temperature [K] (BSCCO) 

FIG. 2. First order phase transition curves of magnetic 
field versus temperature, for YBCO and BSCCO. Parame- 
ters used for YBCO are ai = 2.55, a 2 = 0.01485, <p = 44.1°, 
X ab (0) = 11861 |§, s = 12 A, = ISA, 7 = 5, and 

T c = 92.74 K. Parameters used for BSCCO are ai = 1.0, 
a 2 = 0.00705, cf> = 60°, A ab (0) = 2000A, s = 14A, 
£ a6 (0) = 30A, 7 = 200, and T c = 90 K. For BSCCO we 
use the low field form of the elastic moduli from (hi ) and 
for YBCO we use the high field form. For f we use (2h for 
BSCCO and (g) for YBCO. (For BSCCO we plot B vs. T 
because that is what ref. jn| measured.) The experimental 
points for YBCO come from ref. || and those for BSCCO 
come from ref. fllCf. Inset: Typical A/ versus n. 
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FIG. 4. Surface tension between the vortex solid and liquid 
phases along an interface parallel to the c-axis versus melting 
temperature T m . We use the values of field and temperature 
along the melting curve for YBCO and BSCCO. The surface 
tension is measured in units of Kelvin/sa , where s is the 
spacing between layers and a a = yj <fi a /B is the vortex lat- 
tice spacing. The values of the parameters are the same as in 
Figure 1. 
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FIG. 3. (a) and (b): Entropy jump As per vortex per 
layer versus T m at the transition for YBCO and BSCCO. 
The experimental points for YBCO are from |J and those for 
BSCCO are from Q. (c) and (d): Magnetization jump AM 
versus T m at the first order phase transition for YBCO and 
BSCCO. The experimental points for YBCO are from [^| and 
those for BSCCO are from ^] . For the theoretical points the 
values of the parameters are the same as in Figure 1 for all 
the curves. 
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FIG. 5. Reduced viscosity ri/r} of the vortex liquid versus 
melting temperature T m along the melting curve for YBCO 
and BSCCO. 5 = 1. The values of the parameters are the 
same as in Figure 1. 
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FIG. 6. Shear modulus versus frequency for YBCO and 
BSCCO. For YBCO, H = 5 T, T = 83 K and the defect con- 
centration n = 4 %. For BSCCO, B = 200 G, T — 60 K, and 
n = 5 %. The values chosen are close to those on the melting 
curve. The rest of the values of the parameters are the same 
as in Figure 1. 
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